Search for: All records

Molecular clouds (MCs) are active sites of star formation in galaxies, and their formation and evolution are largely affected by stellar feedback. This includes outflows and winds from newly formed stars, radiation from young clusters, and supernova explosions. High-resolution molecular line observations allow for the identification of individual star-forming regions and the study of their integrated properties. Moreover, state-of-the-art simulations are now capable of accurately replicating the evolution of MCs, including all key stellar feedback processes. We present¹³CO(2–1) synthetic observations of the STARFORGE simulations produced using the radiative transfer code RADMC-3D, matching the observational setup of the SEDIGISM survey. From these synthetic observations, we identified the population of MCs using hierarchical clustering and analysed them to provide insights into the interpretation of observed MCs as they evolve. The flux distributions of the post-processed synthetic observations and the properties of the MCs, namely, radius, mass, velocity dispersion, virial parameter, and surface density, are consistent with those of SEDIGISM. Both samples of MCs occupy the same regions in the scaling relation plots; however, the average distributions of MCs at different evolutionary stages do not overlap on the plots. This highlights the reliability of our approach in modelling SEDIGISM and suggests that MCs at different evolutionary stages contribute to the scatter in observed scaling relations. We study the trends in MC properties, morphologies, and fragmentation over time to analyse their physical structure as they form, evolve, and are destroyed. MCs appear as small diffuse cloudlets in early stages, and this is followed by their evolution to filamentary structures before being shaped by stellar feedback into 3D bubbles and getting dispersed. These trends in the observable properties of MCs are consistent with other realisations of simulations and provide strong evidence that clouds exhibit distinct morphologies over the course of their evolution. 

Stars form in dense cores within molecular clouds, and newly formed stars influence their natal environments. How stellar feedback impacts core properties and evolution has been the subject of extensive investigation. We performed a hierarchical clustering (dendrogram) analysis of a STARFORGE (STAR FORmation in Gaseous Environments) simulation, modelling a giant molecular cloud to identify gas overdensities (cores) and study changes in their radius, mass, velocity dispersion, and virial parameter with respect to stellar feedback. We binned these cores on the basis of the fraction of gas affected by protostellar outflows, stellar winds, and supernovae and analysed the property distributions for each feedback bin. We find that cores that experience more feedback influence are smaller. Feedback notably enhances the velocity dispersion and virial parameter of the cores, more so than it reduces their radius. This is also evident in the linewidth–size relation, according to which cores in higher-feedback bins exhibit higher velocities than their similarly sized pristine counterparts. We conclude that stellar feedback mechanisms, which impart momentum to the molecular cloud, simultaneously compress and disperse the dense molecular gas. 

The ALMA-IMF Large Program provides multi-tracer observations of 15 Galactic massive protoclusters at a matched sensitivity and spatial resolution. We focus on the dense gas kinematics of the G353.41 protocluster traced by N₂H⁺(1−0), with a spatial resolution of ~0.02 pc. G353.41, at a distance of ~2kpc, is embedded in a larger-scale (~8 pc) filament and has a mass of ~2.5 × 10³M_⊙within 1.3 × 1.3 pc². We extracted the N₂H⁺(1−0) isolated line component and decomposed it by fitting up to three Gaussian velocity components. This allows us to identify velocity structures that are either muddled or impossible to identify in the traditional position-velocity diagram. We identify multiple velocity gradients on large (~1 pc) and small scales (~0.2pc). We find good agreement between the N₂H⁺velocities and the previously reported DCN core velocities, suggesting that cores are kinematically coupled with the dense gas in which they form. We have measured nine converging “V-shaped” velocity gradients (VGs) (~20 km s⁻¹pc⁻¹) that are well resolved (sizes ~0.1 pc), mostly located in filaments, which are sometimes associated with cores near their point of convergence. We interpret these V-shapes as inflowing gas feeding the regions near cores (the immediate sites of star formation). We estimated the timescales associated with V-shapes as VG⁻¹, and we interpret them as inflow timescales. The average inflow timescale is ~67 kyr, or about twice the free-fall time of cores in the same area (~33 kyr) but substantially shorter than protostar lifetime estimates (~0.5 Myr). We derived mass accretion rates in the range of (0.35–8.77) × 10⁻⁴M_⊙yr⁻¹. This feeding might lead to further filament collapse and the formation of new cores. We suggest that the protocluster is collapsing on large scales, but the velocity signature of collapse is slow compared to pure free-fall. Thus, these data are consistent with a comparatively slow global protocluster contraction under gravity, and faster core formation within, suggesting the formation of multiple generations of stars over the protocluster’s lifetime. 

The physical mechanisms behind the fragmentation of high-mass dense clumps into compact star-forming cores and the properties of these cores are fundamental topics that are heavily investigated in current astrophysical research. The ALMAGAL survey provides the opportunity to study this process at an unprecedented level of detail and statistical significance, featuring high-angular resolution 1.38 mm ALMA observations of 1013 massive dense clumps at various Galactic locations. These clumps cover a wide range of distances (~2–8 kpc), masses (~10²–10⁴M_⊙), surface densities (0.1–10 g cm⁻²), and evolutionary stages (luminosity over mass ratio indicator of ~0.05 <L/M <450L_⊙/M_⊙). Here, we present the catalog of compact sources obtained with theCuTExalgorithm from continuum images of the full ALMAGAL clump sample combining ACA-7 m and 12 m ALMA arrays, reaching a uniform high median spatial resolution of ~1400 au (down to ~800 au). We characterize and discuss the revealed fragmentation properties and the photometric and estimated physical parameters of the core population. The ALMAGAL compact source catalog includes 6348 cores detected in 844 clumps (83% of the total), with a number of cores per clump between 1 and 49 (median of 5). The estimated core diameters are mostly within ~800–3000 au (median of 1700 au). We assigned core temperatures based on theL/Mof the hosting clump, and obtained core masses from 0.002 to 345M_⊙(complete above 0.23 M_⊙), exhibiting a good correlation with the core radii (M ∝ R^2.6). We evaluated the variation in the core mass function (CMF) with evolution as traced by the clumpL/M, finding a clear, robust shift and change in slope among CMFs within subsamples at different stages. This finding suggests that the CMF shape is not constant throughout the star formation process, but rather it builds (and flattens) with evolution, with higher core masses reached at later stages. We found that all cores within a clump grow in mass on average with evolution, while a population of possibly newly formed lower-mass cores is present throughout. The number of cores increases with the core masses, at least until the most massive core reaches ~10M_⊙. More generally, our results favor a clump-fed scenario for high-mass star formation, in which cores form as low-mass seeds, and then gain mass while further fragmentation occurs in the clump. 

Context.Molecular outflows are believed to be a key ingredient in the process of star formation. The molecular outflow associated with DR21 Main in Cygnus-X is one of the most extreme molecular outflows in the Milky Way in terms of mass and size. The outflow is suggested to belong to a rare class of explosive outflows formed by the disintegration of protostellar systems. Aims.We aim to explore the morphology, kinematics, and energetics of the DR21 Main outflow, and to compare those properties to confirmed explosive outflows in order to unravel the underlying driving mechanism behind DR21. Methods.We studied line and continuum emission at a wavelength of 3.6 mm with IRAM 30 m and NOEMA telescopes as part of the Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE) program. The spectra include (J= 1−0) transitions of HCO⁺, HCN, HNC, N₂H⁺, H₂CO, and CCH, which trace different temperature and density regimes of the outflowing gas at high velocity resolution (~0.8 km s⁻¹). The map encompasses the entire DR21 Main outflow and covers all spatial scales down to a resolution of 3″ (~0.02 pc). Results.Integrated intensity maps of the HCO⁺emission reveal a strongly collimated bipolar outflow with significant overlap of the blueshifted and redshifted emission. The opening angles of both outflow lobes decrease with velocity, from ~80 to 20° for the velocity range from 5 to 45 km s⁻¹relative to the source velocity. No evidence is found for the presence of elongated, “filament-like” structures expected in explosive outflows. N₂H⁺emission near the western outflow lobe reveals the presence of a dense molecular structure, which appears to be interacting with the DR21 Main outflow. Conclusions.The overall morphology as well as the detailed kinematics of the DR21 Main outflow are more consistent with a typical bipolar outflow than with an explosive counterpart. 

Context.Stars form preferentially in clusters embedded inside massive molecular clouds, many of which contain high-mass stars. Thus, a comprehensive understanding of star formation requires a robust and statistically well-constrained characterization of the formation and early evolution of these high-mass star clusters. To achieve this, we designed the ALMAGAL Large Program that observed 1017 high-mass star-forming regions distributed throughout the Galaxy, sampling different evolutionary stages and environmental conditions. Aims.In this work, we present the acquisition and processing of the ALMAGAL data. The main goal is to set up a robust pipeline that generates science-ready products, that is, continuum and spectral cubes for each ALMAGAL field, with a good and uniform quality across the whole sample. Methods.ALMAGAL observations were performed with the Atacama Large Millimeter/submillimeter Array (ALMA). Each field was observed in three different telescope arrays, being sensitive to spatial scales ranging from ≈1000 au up to ≈0.1 pc. The spectral setup allows sensitive (≈0.1 mJy beam⁻¹) imaging of the continuum emission at 219 GHz (or 1.38 mm), and it covers multiple molecular spectral lines observed in four different spectral windows that span about ≈4 GHz in frequency coverage. We have designed a Python-based processing workflow to calibrate and image these observational data. This ALMAGAL pipeline includes an improved continuum determination, suited for line-rich sources; an automatic self-calibration process that reduces phase-noise fluctuations and improves the dynamical range by up to a factor ≈5 in about 15% of the fields; and the combination of data from different telescope arrays to produce science-ready, fully combined images. Results.The final products are a set of uniformly generated continuum images and spectral cubes for each ALMAGAL field, including individual-array and combined-array products. The fully combined products have spatial resolutions in the range 800–2000 au, and mass sensitivities in the range 0.02–0.07M_⊙. We also present a first analysis of the spectral line information included in the ALMAGAL setup, and its potential for future scientific studies. As an example, specific spectral lines (e.g., SiO and CH₃CN) at ≈1000 au scales resolve the presence of multiple outflows in clusters and will help us to search for disk candidates around massive protostars. Moreover, the broad frequency bands provide information on the chemical richness of the different cluster members, which can be used to study the chemical evolution during the formation process of star clusters. 

ABSTRACT The evolutionary sequence for high-mass star formation starts with massive starless clumps that go on to form protostellar, young stellar objects and then compact H ii regions. While there are many examples of the three later stages, the very early stages have proved to be elusive. We follow-up a sample of 110 mid-infrared dark clumps selected from the ATLASGAL catalogue with the IRAM telescope in an effort to identify a robust sample of massive starless clumps. We have used the HCO+ and HNC (1-0) transitions to identify clumps associated with infall motion and the SiO (2-1) transition to identity outflow candidates. We have found blue asymmetric line profile in 65 per cent of the sample, and have measured the infall velocities and mass infall rates (0.6–36 × 10−3 M⊙ yr−1) for 33 of these clumps. We find a trend for the mass infall rate decreasing with an increase of bolometric luminosity to clump mass, i.e. star formation within the clumps evolves. Using the SiO 2-1 line, we have identified good outflow candidates. Combining the infall and outflow tracers reveals that 67 per cent of quiescent clumps are already undergoing gravitational collapse or are associated with star formation; these clumps provide us with our best opportunity to determine the initial conditions and study the earliest stages of massive star formation. Finally, we provide an overview of a systematic high-resolution ALMA study of quiescent clumps selected that allows us to develop a detailed understanding of earliest stages and their subsequent evolution. 

Context.The star formation process leads to an increased chemical complexity in the interstellar medium. Sites associated with high-mass star and cluster formation exhibit a so-called hot core phase, characterized by high temperatures and column densities of complex organic molecules. Aims.We aim to systematically search for and identify a sample of hot cores toward the 15 Galactic protoclusters of the ALMA-IMF Large Program and investigate their statistical properties. Methods.We built a comprehensive census of hot core candidates toward the ALMA-IMF protoclusters based on the detection of two CH₃OCHO emission lines at 216.1 GHz. We used the source extraction algorithm GExt2D to identify peaks of methyl formate (CH₃OCHO) emission, a complex species commonly observed toward sites of star formation. We performed a cross-matching with the catalog of thermal dust continuum sources from the ALMA-IMF 1.3 mm continuum data to infer their physical properties. Results.We built a catalog of 76 hot core candidates with masses ranging from ~0.2M_⊙to ~80M_⊙, of which 56 are new detections. A large majority of these objects, identified from methyl formate emission, are compact and rather circular, with deconvolved full width at half maximum (FWHM) sizes of ~2300 au on average. The central sources of two target fields show more extended, but still rather circular, methyl formate emission with deconvolved FWHM sizes of ~6700 au and 13 400 au. About 30% of our sample of methyl formate sources have core masses above 8M_⊙and range in size from ~1000 au to 13 400 au, which is in line with measurements of archetypical hot cores. The origin of the CH₃OCHO emission toward the lower-mass cores may be explained as a mixture of contributions from shocks or may correspond to objects in a more evolved state (i.e., beyond the hot core stage). We find that the fraction of hot core candidates increases with the core mass, suggesting that the brightest dust cores are all in the hot core phase. Conclusions.Our results suggest that most of these compact methyl formate sources are readily explained by simple symmetric models, while collective effects from radiative heating and shocks from compact protoclusters are needed to explain the observed extended CH₃OCHO emission. The large fraction of hot core candidates toward the most massive cores suggests that they rapidly enter the hot core phase and that feedback effects from the forming protostar(s) impact their environment on short timescales. 

Context. A large fraction of stars form in clusters containing high-mass stars, which subsequently influences the local and galaxy-wide environment. Aims. Fundamental questions about the physics responsible for fragmenting molecular parsec-scale clumps into cores of a few thousand astronomical units (au) are still open, that only a statistically significant investigation with ALMA is able to address; for instance: the identification of the dominant agents that determine the core demographics, mass, and spatial distribution as a function of the physical properties of the hosting clumps, their evolutionary stage and the different Galactic environments in which they reside. The extent to which fragmentation is driven by clumps dynamics or mass transport in filaments also remains elusive. Methods. With the ALMAGAL project, we observed the 1.38 mm continuum and lines toward more than 1000 dense clumps in our Galaxy, withM≥ 500 M_⊙, Σ ≥ 0.1 g cm⁻²andd≤ 7.5 kiloparsec (kpc). Two different combinations of ALMA Compact Array (ACA) and 12-m array setups were used to deliver a minimum resolution of ∼1000 au over the entire sample distance range. The sample covers all evolutionary stages from infrared dark clouds (IRDCs) to H IIregions from the tip of the Galactic bar to the outskirts of the Galaxy. With a continuum sensitivity of 0.1 mJy, ALMAGAL enables a complete study of the clump-to-core fragmentation process down toM∼ 0.3 M_⊙across the Galaxy. The spectral setup includes several molecular lines to trace the multiscale physics and dynamics of gas, notably CH₃CN, H₂CO, SiO, CH₃OH, DCN, HC₃N, and SO, among others. Results. We present an initial overview of the observations and the early science product and results produced in the ALMAGAL Consortium, with a first characterization of the morphological properties of the continuum emission detected above 5σin our fields. We used “perimeter-versus-area” and convex hull-versus-area metrics to classify the different morphologies. We find that more extended and morphologically complex (significantly departing from circular or generally convex) shapes are found toward clumps that are relatively more evolved and have higher surface densities. Conclusions. ALMAGAL is poised to serve as a game-changer for a number of specific issues in star formation: clump-to-core fragmentation processes, demographics of cores, core and clump gas chemistry and dynamics, infall and outflow dynamics, and disk detections. Many of these issues will be covered in the first generation of papers that closely follow on the present publication. 

Context. The electron density ( n e − ) plays an important role in setting the chemistry and physics of the interstellar medium. However, measurements of n e − in neutral clouds have been directly obtained only toward a few lines of sight or they rely on indirect determinations. Aims. We use carbon radio recombination lines and the far-infrared lines of C + to directly measure n e − and the gas temperature in the envelope of the integral shaped filament (ISF) in the Orion A molecular cloud. Methods. We observed the C102 α (6109.901 MHz) and C109 α (5011.420 MHz) carbon radio recombination lines (CRRLs) using the Effelsberg 100 m telescope at ≈2′ resolution toward five positions in OMC-2 and OMC-3. Since the CRRLs have similar line properties, we averaged them to increase the signal-to-noise ratio of the spectra. We compared the intensities of the averaged CRRLs, and the 158 μm-[CII] and [ 13 CII] lines to the predictions of a homogeneous model for the C + /C interface in the envelope of a molecular cloud and from this comparison we determined the electron density, temperature and C + column density of the gas. Results. We detect the CRRLs toward four positions, where their velocity ( v LSR  ≈ 11 km s −1 ) and widths ( σ v  ≈ 1 km s −1 ) confirms that they trace the envelope of the ISF. Toward two positions we detect the CRRLs, and the 158 μm-[CII] and [ 13 CII] lines with a signal-to-noise ratio ≥5, and we find n e −  = 0.65 ± 0.12 cm −3 and 0.95 ± 0.02 cm −3 , which corresponds to a gas density n H  ≈ 5 × 10 3 cm −3 and a thermal pressure of p th  ≈ 4 × 10 5 K cm −3 . We also constrained the ionization fraction in the denser portions of the molecular cloud using the HCN(1–0) and C 2 H(1–0) lines to x (e − ) ≤ 3 × 10 −6 . Conclusions. The derived electron densities and ionization fraction imply that x (e − ) drops by a factor ≥100 between the C + layer and the regions probed by HCN(1–0). This suggests that electron collisional excitation does not play a significant role in setting the excitation of HCN(1–0) toward the region studied, as it is responsible for only ≈10% of the observed emission.